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One of the fascinating properties of the new families of two-dimensional crystals is their high 
stretchability and the possibility to use external strain to manipulate, in a controlled manner, their 
optical and electronic properties. Strain engineering, understood as the field that study how the 
physical properties of materials can be tuned by controlling the elastic strain fields applied to it, has 
a perfect platform for its implementation in the atomically thin semiconducting materials. The ob¬ 
ject of this review is to give an overview of the recent progress to control the optical and electronics 
properties of 2 D crystals, by means of strain engineering. We will concentrate on semiconducting 
layered materials, with especial emphasis in transition metal dichalcogenides (M0S2, WS2, MoSe2 
and WSe2). The effect of strain in other atomically thin materials like black phosphorus, silicene, 
etc., is also considered. The benefits of strain engineering in 2 D crystals for applications in nano¬ 
electronics and optoelectronics will be revised, and the open problems in the field will be discussed. 
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I. INTRODUCTION 

The isolation of graphene by mechanical exfoliation 
in 2004 [T] opened the door to study a broad fam¬ 
ily of two-dimensional materials ms, almost unex¬ 
plored at that time. These materials are character¬ 
ized by strong covalent in-plane bonds and weak in¬ 
terlayer van der Waals interactions which give them 
a layered structure. This rather weak interlayer in¬ 
teraction can be exploited to extract atomically thin 
layers by mechanical [2] or chemical exfoliation meth¬ 
ods [6]. In the last five years a wide variety of ma- 
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terials with very different electronic properties (rang¬ 
ing from wide-bandgap insulators to superconductors) 
have been explored [MU. These materials are ex¬ 
pected to complement graphene in applications for 
which graphene does not possess the optimal proper¬ 
ties. Two-dimensional semiconductors with an intrin¬ 
sic large bandgap are a good example of materials that 
could complement graphene EE]. While the out¬ 
standing carrier mobility of graphene makes it very at¬ 
tractive for certain electronic applications (e.g. high- 
frequency electronics), the lack of a bandgap dramat¬ 
ically hampers its applicability in digital electronics. 
Graphene field-effect transistors suffer from large off- 
state currents due to the low resistance values at the 
charge neutrality point. The large off-state current 
also translates to a large dark current in graphene- 
based photodetectors which limits their performance. 
The handicaps of gapless graphene-based electronic 
devices have motivated experimental efforts towards 
opening a bandgap in graphene in order to combine 
its excellent carrier mobility with a sizable bandgap. 
Several methods have been employed: lateral con¬ 
finement [HHni, application of a perpendicular elec¬ 
tric field in bilayer graphene nans], hydrogenation 
[2Qf[23] . or controlled nanoperforation [24f[26] among 
others. However, in all cases the magnitude of the 
bandgap is not enough to ensure high-performance 
room-temperature and/or the mobility of the treated 
graphene severely decreases. 

The difficulty of achieving high-mobility gapped 
graphene has triggered the study of semiconducting 
analogues to graphene in materials with an intrinsic 
bandgap [SIIITHSI]. For instance single-layer M 0 S 2 , 
the most studied 2D semiconductor so far, presents a 
direct bandgap (1.8 eV) and a large in-plane mobil¬ 
ity (up to 200 cm^V“^s“^) [32H37] . The intrinsic 

bandgap of M 0 S 2 is essential for many applications, 
including transistors for digital electronics or certain 
optoelectronic applications. For instance, field-effect 
transistors based on single-layer M 0 S 2 showed room- 
temperature current on/off ratios of 10^ — 10^ and 
ultralow standby power dissipation [14]. Logic cir¬ 
cuits, amplifiers and photodetectors based on mono- 
layer M 0 S 2 have also been recently demonstrated 

jsSlEHlEl]. 

Apart from their outstanding figures-of-merit and 
electrical performances, atomically thin semiconduc¬ 
tors have also shown very interesting mechanical prop¬ 
erties m, unmatched by conventional 3D semicon¬ 
ductors mnn. Nanoindentation experiments with 
an atomic force microscope tip on freely suspended 
single-layer M 0 S 2 have shown mechanical properties 


approaching those predicted by Griffith for ideal brit¬ 
tle materials in which the fracture point is dominated 
by the intrinsic strength of their atomic bonds and 
not by the presence of defects m- While the ideal 
breaking stress value is expected to be one ninth of 
the Youngs modulus, for single layer M 0 S 2 h reaches 
^ 1/8 im |44]. This almost ideal behavior is at¬ 
tributed to a low density of defects on the fabricated 
devices, probably due to the lack of dangling bonds 
or other surface defects and their high crystallinity. 
This is in clear contrast with conventional 3D semi¬ 
conductors. Indeed, while silicon typically breaks at 
strain levels of ^ 1.5%, M 0 S 2 does not break until 
> 10% strain levels and it can be folded and wrin¬ 
kled almost at will |45j [46]. On the other hand, the 
elastic properties of 2D crystals at long wavelengths 
are dominated by anharmonic out-of-plane vibrations, 
and a controlled creation of vacancies [47] or strain 
[48l l49] can be used to enhance the strength of these 
materials. Finally, the ultimate crystal instability (on¬ 
set of negative phonon frequencies) of several families 
of two-dimensional materials was investigated in Ref. 

m- 

This outstanding stretchability of 2D crystals 
promises to revolutionize the field of strain engineer¬ 
ing due to the unprecedented tunability levels pre¬ 
dicted for these materials that could lead to ”strain- 
tronic” devices devices with electronic properties that 
are engineered through the introduction of mechan¬ 
ical deformations. Single-layer M 0 S 2 , for instance, 
is expected to undergo a direct-to-indirect bandgap 
transition at ^ 2% of tensile uniaxial strain, and a 
semiconducting-to-metal transition at 10 — 15% of ten¬ 
sile biaxial strain [5TH53] . Such huge sensitivity of 
the bandgap (from 1.8 eV to 0 eV) can be realisti¬ 
cally exploited in M 0 S 2 and other transition metal 
dichalcogenides, since these materials can withstand 
such levels of strain without rupture [45l [46]. This 
behavior has to be compared with the poor window 
of tunability of only 0.25 eV achieved for strained sil¬ 
icon at 1.5% biaxial strain, due to the breaking of 
the bonds beyond that value of strain [54]. Another 
game-changing feature of atomically thin semiconduc¬ 
tors is that one can deform them only at a specific lo¬ 
cation making it possible to generate localized strain 
profiles. One possible advantage of local strain engi¬ 
neering has been recently proposed by Feng et al. in 
Ref. [52] when a M 0 S 2 membrane is subjected to a 
point deformation at its centre, a continuous bend¬ 
ing of the energy levels of electrons, holes and ex- 
citons is produced. Upon illumination electron-hole 
pairs are generated and the semiclassical potential for 
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the quasiparticles exerts forces on the electron and 
holes in opposite directions. However, the large exci- 
ton binding energy will force the electron-hole pairs to 
hold together and migrate towards the centre of the 
membrane. Therefore, localized strains can be used to 
concentrate excitons in a small region of the semicon¬ 
ductor crystal with interest for fundamental physics 
(e.g. exciton condensates) and applications (e.g. so¬ 
lar energy funnel). 

The present review covers the different levels in the 
study of strain engineering in 2D crystals, starting 
from the fabrication in Sec. EAl The different tech¬ 
niques to induce strain are discussed in Sec. IIB and 
their effects on the optical and electronic properties 
are summarized in Sec. |II C| where we list the main 
experimental results. From a theoretical point of view, 
the electronic structure of semiconducting 2D mate¬ 
rials is revised in Sec. [mil In Sec. jm^ we go 
over the effect of homogeneous strain on these ma¬ 
terials, focusing on transition metal dichalcogenides, 
black phosphorus and silicene, and the effect of non- 
uniform strain is discussed in Sec. Imcl In Sec. lEI 
we give our conclusions and final remarks. 


II. TWO DIMENSIONAL CRYSTALS 





A. Fabrication 

In this section the main fabrication methods to 
isolate or synthesize 2D materials will be briefly de¬ 
scribed. For a more detailed discussion about fabrica¬ 
tion methods we address the reader to recently pub¬ 
lished review articles that are focused on this specific 
topic isa ES]. Different fabrication techniques typi¬ 
cally yield different size, thickness and even quality of 
the fabricated 2D materials and therefore the selec¬ 
tion of the isolation method strongly depends on the 
desired application. 


FIG. 1 : Isolation of single-layer M0S2 by differ¬ 
ent techniques, (a) Optical image of a M0S2 flake 
deposited on a Si02/Si substrate by mechanical exfoli¬ 
ation. (b) Optical image of a single-layer M0S2 crystal 
grown on a Si 02 /Si by chemical vapour deposition method 
(adapted from Ref. [ 65 ] with permission of Nature Pub¬ 
lishing Group), (c) Transmission electron microscopy im¬ 
age of a M0S2 flake isolated by ultrasonication of M0S2 
powder in an organic solvent. The inset shows a high- 
resolution scanning transmission electron microscopy im¬ 
age of the crystal lattice of the single layer region. Adapted 
from Ref. [ 73 ] with permission of Science. 


1. Mechanical exfoliation 

The weak van der Waals interaction between the 
layers makes it possible to cleave thin crystalline flakes 
by peeling off the surface of a bulk layered material 
that is adhered to a piece of sticky tape. This method 
to isolate thin sheets from a bulk layered material is re¬ 
ferred to in the literature as Scotch tape method, me¬ 
chanical exfoliation or micromechanical cleavage and 
it has proven to be a simple yet powerful technique to 


obtain high-quality two-dimensional sheets. The pro¬ 
cess is based on repeatedly cleaving a bulk layered ma¬ 
terial with a sticky tape to generate a large quantity 
of thin crystallites [1] . These crystallites can be trans¬ 
ferred to an arbitrary flat substrate by gently pressing 
the tape against the surface of the substrate and peel¬ 
ing it off slowly. More details on this technique can 
be found in the pioneering work of Novoselov, Geim 
et al [ 2 ]. 

The main limitation of this technique is the fact 
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(C) 

A 2D material is deposited on the 
center of a elastic substrate 



FIG. 2 : Uniaxial straining of 2 D materials with bending and tensile test geometry, (a) Cartoon of a flexible 
substrate with a 2 D material deposited onto its surface. By using two-point bending test geometry, the substrate can 
be bent to apply a controlled uniaxial tensile strain to the 2 D material, (b) Simplified scheme of the cross section of the 
flexible substrate before and after bending it. (c) Scheme of the approach to apply uniaxial tensile strain by elongating 
an elastic substrate. Panel (a) adapted from [ 80 ] with permission of the American Physical Society. 


that during the deposition step flakes with various 
thicknesses are transferred all over the surface and 
only a small fraction of them are few-layers thick. 
Nonetheless, this issue can be overcame to a large ex¬ 
tent because optical microscopy can be used to iden¬ 
tify thin 2D flakes in a fast and reliable way m\- Fig. 
I^a) shows an example of optical microscopy image 
of a single-layer M 0 S 2 Hake. Its thickness can be dis¬ 
tinguished at glance by its optical contrast, due to 
a combination of interference color and optical ab¬ 
sorption. Several quantitative studies of the optical 
contrast of many 2D materials, providing a guide to 
optically identify 2D materials, have been recently re¬ 
ported: M 0 S 2 laiSimO], NbSe2,[9l|60], WSe 2 EHEOj, 
TaS 2 [SS], TaSe 2 [61], mica [m[ 62 ] . 


2. Chemical vapour deposition or epitaxial growth 

2D materials can be also prepared by bottom-up 
approaches where layers are synthesized by assem¬ 
bling their elemental constituents. Transition metal 
dichalcogenides (TMDs) have been synthesized by 
vapour transport method {63] [64], chemical vapour 
deposition [65H68] or van der Waals epitaxy [69] - 
lll]. High-quality molybdenum and tungsten dichalco¬ 
genides can be directly grown on insulating substrates 
with large single-crystal domains of 10 to 10000 
|6l|65l|68]. Fig. Efb) shows an optical microscopy 
image of a typical single-crystalline M 0 S 2 monolayer 
grown by chemical vapour deposition [65] . 


3. Liquid-phase exfoliation 

Alternatively to the approaches described below, 
atomically thin 2D crystals can be also isolated by ex¬ 
foliating bulk layered materials immersed in a liquid 
medium using two main approaches: direct sonication 
in a solvent or intercalation/expansion/sonication. 
We address the reader to Refs. [MIS] and to Refs. 
imiTE] for details about these two liquid phase exfo¬ 
liations approaches, respectively. 

Due to the combination of high yield and low cost, 
this technique is a prospective fabrication approach to 
fabricate large quantities of atomically thin crystals 
but it yields small size flakes (0.2 — and it lacks 

of sufficient control to achieve monodisperse solutions. 
Note that liquid phase exfoliated flakes have not been 
used in the strain engineering studies reported so far, 
probably because of their reduced lateral sizes. Fig. 
[^c) shows a transmission electron microscopy image 
of a M0S2 flake fabricated by liquid phase exfoliation 
of M0S2 powder [73] . 


B. Techniques to induce strain 

In this section the methods employed to strain en¬ 
gineer 2D materials will be described. Leaving aside 
techniques which have been employed for specific pur¬ 
poses in graphenes and not yet applied for the fami¬ 
lies of 2D compounds we focus here (transition metal 
dichalcogenides, black phosphorus, silicenes, ...), we 
will address in the following four techniques that have 
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already proven to be feasible and currently employed 
for these materials: bending of flexible substrates, 
elongating an elastic substrate, piezoelectric stretch¬ 
ing, substrate thermal expansion and controlled wrin¬ 
kling. These methods will be described in separated 
subsections. A summary of this section will com¬ 
pare the results obtained by the different straining 
approaches. 


1. Bending of a flexible substrate 

Controllable uniaxial tensile strain can easily be ap¬ 
plied to 2D materials by using a modified version of 
a bending test. The 2D material is deposited onto a 
flexible substrate (directly by mechanical exfoliation 
of a bulk layered material or by transferring a vapour- 
transport growth film). By bending the flexible sub¬ 
strate as depicted in Fig. [^a) the topmost surface is 
stretched, transferring a uniaxial tensile strain to the 
2D material that lays on the surface. The bending 
of the flexible substrate is typically carried out with 
a micrometer stage in a two-points, three-points or 
four-points bending geometry. 

The strain achieved by bending the flexible sub¬ 
strate can be obtained from a continuum mechanics 
model for elastic beams assuming that the radius of 
curvature (R) due to bending the substrate is larger 
than the substrate thickness (h) [79] (see sketch in 
Fig. [^b)). Also, the dominant deformation of the 
beam should be in the longitudinal direction. This 
means that shear stresses and stresses normal to the 
neutral axis are negligible. Based on these assump¬ 
tions the strain can be obtained as m- 

s = h/2R (1) 

The van der Waals interaction between the 2D ma¬ 
terial and the surface of the flexible substrate is typi¬ 
cally strong enough to clamp the nanosheet and pre¬ 
vent slippage for strain levels up to 1%. For higher 
strain levels, the sheets tend to slip and the tension 
is released hampering the acquisition of reproducible 
datasets. In order to reduce the slippage issue, one 
can evaporate metal strips onto the nanosheets to act 
as clamping points [8T] . 

This technique to apply strain has been successfully 
employed to study the changes in the Raman vibra¬ 
tional modes of graphene subjected to tensile strain 
[sol [82]-[84]. More recently, this straining technique 
has been adopted to study the role of uniaxial de¬ 
formation in the optical properties of atomically thin 





FIG. 3 : Uniform biaxial straining by using piezo¬ 
electric substrates, (a) Schematic diagram of the me¬ 
chanical deformation of a piezoelectric material when an 
electric field is applied, (b) Cartoon of the experimen¬ 
tal setup employed in Ref. [ 89 ] to apply uniform biaxial 
strain to trilayer M0S2 (adapted with permission of the 
American Chemical Society). 


M 0 S 2 which will be described in the next section [85]- 
[87] . The maximum reported uniaxial strain achieved 
by this technique ranges from 0.5 % to 2.5%, depend¬ 
ing on the exact experimental conditions. 


2. Elongating the substrate 

Alternatively to the previous technique homoge¬ 
neous uniaxial tensile strain can be also applied by 
depositing the 2D material on an elastic substrate and 
elongating it with a straining stage, like the ones used 
for tensile tests in material science (see Fig. ic)). 
This method has been recently used to apply uniaxial 
tensile strain (up to 4.0%) to monolayer WS 2 [88] . 


3. Piezoeleetrie stretehing 

Piezoelectric materials are materials that respond 
to external electric fields with a mechanical defor¬ 
mation. Fig. I^a) shows a schematic diagram of a 
piezoelectric material that is elongated in the verti- 
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cal direction and shrink in the horizontal one when 
an electric field is applied (by applying a bias volt¬ 
age between two faces of the piezoelectric material). 
Therefore, the use of piezoelectric substrates is very 
appealing to apply controllable strains to 2D materi¬ 
als as one can easily regulate the deformation of the 
substrate by simply varying the applied voltage to the 
piezoelectric material. 

Hui et al have transferred trilayer M 0 S 2 onto a 
[Pb(Mgi/ 3 Nb 2 / 3 ) 03 ]o. 7 -[PbTi 03 ]o .3 piezoelectric sub¬ 
strate [89] . By applying a bias voltage to the substrate 
(as illustrated in the sketch in Fig. |^b)) the surface 
where the flake is deposited shrink, subjecting the tri¬ 
layer M 0 S 2 to a uniform compressive biaxial strain. 
With this technique only moderate compressive biax¬ 
ial strain levels have been achieved (a maximum 0.2 
% have been reported so far). Note, however, that 
biaxial strain has a stronger effect on the band struc¬ 
ture of transition metal dichalc ogenide materials than 
uniaxial strain (see Sec. IIC 2). 


4- Exploiting the thermal expansion mismateh 

Thermal expansion mismatch is typically seen as a 
drawback in microfabrication of thin film devices pro¬ 
cess as it leads to undesired strain during the different 
thermal processes needed for the fabrication of micro¬ 
devices. One can, however, exploit the thermal expan¬ 
sion mismatch to engineer strain on a semiconducting 
material in order to modify its optoelectronic prop¬ 
erties in a controllable way. The sketch in Fig. [^a) 
shows an example of how a large mismatch between 
the thermal expansion of a 2D semiconducting mate¬ 
rial and of the substrate where it is laying can be ex¬ 
ploited to induce uniform biaxial tensile strain. While 
2D materials typically have negative thermal expan¬ 
sion coefficients, the substrate can be intentionally se¬ 
lected to have a very large positive thermal expansion 
coefficient. Therefore, by warming up the substrate 
the 2D material will experience a tensile biaxial strain. 

This effect has been observed by Plechinger et al. 
in single-layer M 0 S 2 sandwiched between two dielec¬ 
tric layers m- By changing the sample temperature 
a shift in the photoluminescence spectra evidenced a 
change in the electronic band structure that can be 
explained by the tension induced during the thermal 
cycling by the mismatch in thermal expansion coeffi¬ 
cients. In these experiments, however, the magnitude 
of the induced strain was modest because of the rela¬ 
tively small thermal expansion coefficient of the dielec¬ 
tric layers employed. More recently, it has been shown 
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Laser-induced 
local heating 


FIG. 4: Uniform biaxial straining by exploiting the 
thermal expansion of the substrate, (a) Schematic 
diagram of the process to strain atomically thin materials 
by using substrates with strong adhesion and very large 
thermal expansion coefficient, (b) Cartoon of the strain¬ 
ing process using a focalized laser as a heating source. 
Adapted from Ref. m with permission of the Institute 
Of Physics. 


how this effect can be magnified by selecting a sub¬ 
strate with a large thermal expansion coefficient. For 
example, by using poly-dimethil siloxane instead of 
conventional dielectric substrates it has been recently 
achieved 0.23 % of biaxial tensile strain by increas¬ 
ing the temperature of the sample by 150 °C m- 
In that reference, they also demonstrated that this 
technique can be modified to attain non-homogeneous 
strain profiles. In fact Plechinger et al. proposed 
to employ a focalized laser beams to locally increase 
the temperature of the substrate at certain positions 
which would lead to more complex straining patterns 

m- 


5. Controlled wrinkling 

Transferring 2D materials onto elastomeric sub¬ 
strates allows one to generate non-homogeneous de¬ 
formations in a relatively simple way by exploiting 
the phenomena known as buckling-induced rippling 
or wrinkling. Fig. [^a) shows a schematic diagram 
of the buckling-induced rippling/wrinkling process. A 
thin film is deposited onto an elastomeric substrate 
that has been uniaxially deformed prior to the deposi¬ 
tion step. Once the thin film is deposited the flexible 
substrate strain is released which lead to a uniaxial 
compression of the thin film. If the initial strain is 
below a critical value, the thin film will deform upon 
the compressive strain in a sinusoidal way without de¬ 
laminating from the substrate (rippling). For higher 
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FIG. 5: Non-homogeneous straining by wrinkling 2D materials, (a) Schematic diagram of the fabrication 
process of wrinkled nanolayers. An elastomeric substrate is stretched prior depositing the 2D material by mechanical 
exfoliation. The strain is then released producing buckling-induced delamination of the flakes, (b) and (c) High angle 
scanning electron microscopy images of wrinkled M 0 S 2 and ReSe 2 , respectively. Panel (a) adapted from Ref. |94] with 
permission of the American Chemical Society. Panels (b) and (c) are adapted from Ref. [93] and Ref. [95] respectively 
with permission of the American Chemical Society. 


initial strain levels, the thin film delaminates from the 
substrate forming wrinkles that accumulate the de¬ 
formation, separated by fiat unstrained regions. This 
process has been previously used to generate periodic 
ripples in graphene [saES] and it has been recently 
used to generate non-homogeneous strain profiles in 
atomically thin M 0 S 2 and ReSe 2 [MlllS]- Fig. |^b) 
and (c) show high angle scanning electron microscopy 
images of wrinkles generated in thin M 0 S 2 and ReSe 2 
layers using the buckling-induced delamination pro¬ 
cess. This method provides a powerful, yet simple, 
way to achieve strongly localized strains up to a 2.5% 
[94] . As it will be discussed in the coming sections, 
localized strains can be exploited to collect photogen¬ 
erated charge carriers by generating an exciton funnel. 

The maximum local strain e is induced on the top¬ 
most part of the winkles and it can be estimated by 
the following formula 5 ^ 7r^M/(l — z>^)A^ [94] [96] . 
where v is the Poissons ratio, h is the thickness of the 
flake, and S and A are the height and width of the 
wrinkle extracted from atomic force microscopy. Al¬ 
ternatively to the strain estimation from the geometry 
of the wrinkles, the maximum strain level can also be 
obtained from the shift of the Raman peak. In 
fact. Rice et al. studied the effect of uniaxial tensile 
strain on single-layer M0S2, finding that the Ra¬ 


man peak is shifted towards lower Raman shifts with 
uniaxial tensile strain (it shifts by — 1.7cm“^ per %) 
while the Ai^ peak is almost strain independent m- 

6. Comparison between different teehniques to induee 
strain 

In this section, we summarize the different straining 
techniques described in the previous sections in order 
to facilitate a comparison by the reader. This sum¬ 
mary allows one to find the most adequate straining 
technique for each application. Table |T| displays use¬ 
ful information for each straining method such as the 
type of strain achieved (uniaxial/biaxial and homo¬ 
geneous/inhomogeneous), the maximum strain level 
reported for each technique in the literature and the 
2D materials that have been strain engineered using 
these methods. 


C. Bandgap engineering 

In this section the recent experimental progress on 
modifying the band structure of atomically thin semi¬ 
conducting materials, beyond graphene, by strain en- 









Straining technique 

Type of strain 

Max. strain 

Material 

Ref. 




2 . 4 % 

M 0 S 2 

185J 

Bending of a flexible 
substrate 

Uniaxial 

Homogeneous 

0 . 5 % 

0 . 8 % 

M 0 S 2 

M 0 S 2 

M 

m 




2 . 1 % 

WSe2 

EH] 

Elongating the substrate 

Uniaxial 

Homogeneous 

4 . 0 % 

WS 2 

[88] 

Piezoelectric stretching 

Biaxial 

Homogeneous 

0 . 2 % 

M 0 S 2 

[89] 

Exploiting the thermal ex¬ 

pansion mismatch 

Biaxial 

Homogeneous 

0.23% 

M 0 S 2 

[91] 

Controlled wrinkling 

Uniaxial 

Inhomogeneous 

2.5% 

1.6% 

M 0 S 2 

ReSe2 

194] 

m 


TABLE I: Comparison between the different techniques to apply strain. The corresponding bibliographic reference is 
given in square brackets in the last column. 


gineering will be reviewed. The first experimental 
works that demonstrate the control over the optoelec¬ 
tronic properties of atomically thin semiconducting 
materials are relatively young and are almost limited 
to M 0 S 2 . In 2012, Plechinger et al. reported a shift 
of the photoluminescence yield of single-layer M 0 S 2 
sandwiched between two dielectric layers when the 
temperature of the stack was cycled between 4 K and 
240 K [90]. Similar results were found by Yan et al. 
[99] for single-layer M 0 S 2 samples also sandwiched be¬ 
tween two dielectric layers and heated with an increas¬ 
ingly high power laser. These results were understood 
in the basis that the thermal expansion coefficient mis¬ 
match between the M 0 S 2 and the surrounding dielec¬ 
tric layer generated a biaxial strain when the tempera¬ 
ture of the stack was varied. In these works, however, 
the magnitude of the applied strain was modest and 
the maximum attained strain level was not quantified. 
Soon after these works, five different groups simul¬ 
taneously reported their experimental result demon¬ 
strating the control over the band gap of atomically 
thin M 0 S 2 by using mechanical strain [85H871 [89l l94] . 
These works employed complementary straining tech¬ 
niques and provided a quantitative estimation of the 
strain level applied: bending of a flexible substrate 
I85H87] . using of a piezoelectric substrate [89], and 
controllably wrinkling [94|. More recently, the use of 
substrates with very large thermal expansion coeffi¬ 
cient have been used to engineer biaxial tensile strain 
in single-layer M 0 S 2 m- We refer the reader to pre¬ 
vious section to find a detailed discussion about these 
different techniques to apply strain. 


1. Homogeneous uniaxial strain 

He et al.^ Conley et al. and Zhu et al. transferred 
single- and bilayer M 0 S 2 onto flexible substrates and 
using a bending test experimental configuration they 
subjected these layers to homogeneous tensile uniaxial 
strain [85ll87] . In these works it was found that the 
excitonic features, present in the photoluminescence 
emission and absorption spectra, linearly shift towards 
lower energy when the tensile strain is increased, in¬ 
dicating a reduction in the band gap upon homoge¬ 
neous tensile uniaxial strain, as shown in Fig. The 
magnitude of the photoluminescence and/or absorp¬ 
tion shifts upon strain reported in Refs. [85ll87] are in 
relatively good agreement and discrepancies might be 
attributed to a difference in the Poissons ratio of the 
different flexible substrates employed in the different 
works, as the bending of a beam induces not only an 
in-axis tension but also an off-axis compression whose 
magnitude is dictated by the substrate Poisson’s ratio. 
Uniform uniaxial tensile strain has been also employed 
to tune the band structure of other 2D semiconduc¬ 
tors beyond M 0 S 2 . In Refs. [88l [98] a detailed analysis 
of the photoluminescence emission of WS 2 and WSe 2 
nanosheets upon uniaxial tensile strain is reported. 


2. Homogeneous biaxial strain 

Hui et al. employed a piezoelectric substrate 
to apply homogeneous compressive strain to a tri¬ 
layer M0S2 flake, mechanically exfoliated onto a 
[Pb(Mgi/ 3 Nb 2 / 3 ) 03 ]o. 7 -[PbTi 03 ]o ,3 substrate [89]. A 
CVD grown graphene layer was transferred on top of 
the trilayer M0S2 flake as transparent and flexible top 
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FIG. 6: Band engineering with homogeneous uni¬ 
axial strain, (a) Photoluminescence spectra of a mono- 
layer M 0 S 2 as it is strained from 0 to 1.8%. (inset) 
Schematic of the beam bending apparatus used to strain 
M 0 S 2 . (b) Strain dependent position of the A peak for 
several monolayers (colors represent different samples). 
Adapted from Ref. [85] with permission of the American 
Chemical Society. 


electrode to apply a uniform bias voltage to the piezo¬ 
electric material. The changes in the electronic band 
structure, induced by the straining with the piezo¬ 
electric substrate, were also monitored by the shift 
of the A excitonic feature in the photoluminescence 
spectra. Fig. [^a) shows a collection of photolumines¬ 
cence spectra acquired at different compressive strain 
levels. The position of the photoluminescence peak, 
due to the generation of the A exciton, as a function 
of the strain is displayed in Fig. [^b). From this mea¬ 
surement Hui et al determined a tunability of the A 
exciton of trilayer M 0 S 2 of +300 meV/% of compres¬ 
sive biaxial strain. 

As discussed in the previous section, another 
method to apply homogeneous biaxial strain to 2D 


(a) 



(b) Energy (eV) 



FIG. 7: Band engineering with homogeneous bi¬ 
axial strain, (a) Photoluminescence spectra of trilayer 
M 0 S 2 as it is strained with a piezoelectric substrate, (b) 
Strain dependent position of the A peak. Adapted from 
Ref. [89] with permission of the American Ghemical Soci¬ 
ety. 


semiconductors is to transfer them to a substrate with 
a large thermal expansion coefficient and heating up 
the substrate. In this way the substrate will expand 
at a higher rate than the 2D semiconductor material 
and thus (if the adhesion to the substrate is strong 
enough to avoid slippage) a net biaxial tension will 
be applied to the 2D material. This straining tech¬ 
nique has been employed by Plechinger et al to strain 
monolayer M 0 S 2 . By using a poly-dimethil siloxane 
substrate instead of conventional dielectric substrates 
large biaxial strain can be achieved m- Note that the 
thermal expansion coefficient of poly-dimethil siloxane 
is 1000 times larger than that of Si02 and 150 times 
larger than that of bulk M 0 S 2 . By heating up the 
substrate to 150 °C a biaxial tensile strain of 0.23% 
was achieved, yielding a shift of the A exciton of —60 
meV. By measuring on a Si02 substrate (which has a 
thermal expansion smaller than M 0 S 2 and a poor ad¬ 
hesion) as a reference substrate Plechinger et al de- 
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FIG. 8: Band engineering with inhomogeneous 

uniaxial strain, (a) Photoluminescence spectra of a tri¬ 
layer M 0 S 2 acquired at two different locations: flat (0% 
strain) and on top of a wrinkle (1.1% strain), (b) Strain 
dependent position of the A exciton peak. According to a 
tight binding calculation, including a correction to account 
for the exciton funnelling, the expected data lay within the 
light-blue area. Adapted from Ref. m with permission 
of the American Chemical Society. 


termined that 40% of the observed shift is due to the 
substrate induced biaxial tension and the other 60% 
is due to the intrinsic expansion of M 0 S 2 upon heat¬ 
ing. Therefore, a strain tunability of —105 nieV/% of 
tensile biaxial strain was found in these experiments. 
Similar results were also obtained by employing a fo¬ 
calized laser to increase the temperature of the sub¬ 
strate instead of a hotplate and the use of focalized 
laser heating to engineer more complex straining pat¬ 
terns was also explored in Ref. IS]. 
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FIG. 9: Summary of the state-of-the-art on exper¬ 
iments and theoretical predictions on strain engi¬ 
neering in atomically thin M0S2, WSe2 and ReSe2. 

The hgure also includes some theoretical predictions of 
the band gap change in single-layer M 0 S 2 upon very large 
strains for comparison. These works will be discussed in 
more detail in Sec. nni 
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3. Non-homogeneous uniaxial strain 


Alternatively to the previous experimental ap¬ 
proaches to apply strain, that employed uniform 
straining techniques, Castellanos-Gomez et al. em¬ 
ployed a controllable wrinkling method that make it 
possible to achieve highly localized strains [94]. Win¬ 
kles are fabricated in the M 0 S 2 layers by transfer¬ 
ring ultrathin M 0 S 2 flakes onto a pre-stretched poly- 
dimethil siloxane and subsequently releasing the strain 
(see the previous section for more details). 

The photoluminescence emission of the wrinkled 
M 0 S 2 samples depends on the position. On the flat re¬ 
gions of the sample (no strain) the photoluminescence 
shows a prominent peak corresponding to the genera¬ 
tion of the A exciton (see Fig. |^a)). On the topmost 
part of the wrinkles (where the uniaxial strain due to 
the bending of the layer is maximum) the A exciton 
peak shifts towards lower energies as shown in Fig. 
[^a). This local variation of the photoluminescence 
indicates that the band structure spatially varies fol¬ 
lowing the strain profile on the sample. This can be 
exploited to attract excitons, generated far from the 
maximum strain locations, towards the topmost part 
of the wrinkles m- This concept, referred to in the 
literature as funneling effect, has been theoretically 











11 



Work 

Material 

Experiment 

Strain tunability 

Homogeneous 

uniaxial 

Conley et al [85] 

IL M 0 S 2 

PL (A exciton) 

-45 ± 7 meV/% 

2L M 0 S 2 

PL (A exciton) 

PL (I exciton) 

-53 ± 10 meV/% 
-129 ±20 meV/% 

He et al. [86] 

IL M 0 S 2 

PL & Abs (A exciton) 
PL & Abs (B exciton) 

-64 ± 5 meV/% 
-68 ± 5 ineV/% 

2L M 0 S 2 

PL (A exciton) 

Abs (A exciton) 

PL & Abs (B exciton) 
PL (I exciton) 

-48 ± 5 ineV/% 
-71 ± 5 meV/% 
-67 ± 5 meV /% 
-77 ± 5 meV/Vo 

Zhu et al. [87] 

IL M 0 S 2 

PL (A exciton) 

—48 meV/% 

2L M 0 S 2 

PL (A exciton) 

PL (I exciton) 

—46 meV/% 

—86 meV/% 

Wang et al. [88] 

IL WS 2 

PL (A exciton) 

PL (A“ exciton) 

PL (indirect) 

— 11 meV/% 

— 11 meV/% 

— 19 meV/% 

Desai et al. [98] 

1L-4L WSe2 

PL 

N/A 

Homogeneous 

biaxial 

Hui et al. [89] 

Plechinger et al. [91] 

3L M 0 S 2 

PL (A exciton) 

±300 meV/%* 

IL M 0 S 2 

PL (A exciton) 

— 105 meV/% 

Inhomogeneous 

uniaxial 

Castellanos-Gomez et al. [94] 

3L-5L M 0 S 2 

PL (A exciton) 

—60 meV/% 

Yang et al. [9^ 

1L-2L ReSe2 

PL (A exciton) 

-37 meV/% 


TABLE II: Comparison between the different strain engineering experimental results available in the literature for 2D 
semiconducting materials. 


proposed by Feng et al [52] and it will be discussed 
in more detail in the theory sections of this review 
article. 

Fig. [^b) shows the strain dependence of the A- 
exciton shift, measured in several wrinkles with dif¬ 
ferent maximum strain levels. Unlike the experimen¬ 
tal data acquired with the previous approaches, the 
A-exciton shift shows a marked sublinear trend and a 
relatively large dispersion (larger than that expected 
from the measurement setup). A tight binding calcu¬ 
lation shows that these two features can only be re¬ 
produced by including a correction in the theory that 
accounts for the drift of the photogenerated excitons 
towards regions with larger strain before recombina¬ 
tion. Therefore, the presence of these two features can 
be considered the smoking gun of the funneling effect 
taking place in the inhomogeneously strained M 0 S 2 
layer. 

The inhomogeneous uniaxial straining has been also 
applied to modify the bandstructure of ReSe 2 by Yang 
et al. [95|. In this work, the role of strain on the 
electrical and magnetic properties of ReSe 2 is also ex¬ 
plored using the same wrinkling method. 


4 . Summary of the different hand engineering 
experiments 

In this section, we present a summary of the dif¬ 
ferent strain engineering experiments carried out on 
2D semiconducting materials. Table [11] displays the 
strain tunability determined for different materials 
(M 0 S 2 , WS 2 , WSe 2 and ReSe 2 ) employing different 
straining techniques. Table |TT| also displays the kind of 
experimental measurement (photoluminescence spec¬ 
troscopy or absorption spectroscopy) as well as the 
spectroscopic feature studied (A, A“, B or I exciton). 

Fig. 1^ shows a graphic summary of the state-of- 
the-art on experiments on strain engineering in atom¬ 
ically thin 2D semiconductors. The Figure displays 
the maximum band gap change (measured with op¬ 
tical spectroscopy techniques) versus the maximum 
strain load achieved in each experimental work. Note 
that, in Fig. [^experimental results obtained for dif¬ 
ferent materials, with different straining techniques 
and applying different kind of strains are all plotted 
together. Some theoretical predictions of the band 
gap change for larger strain loads are also included 
in the Figure to illustrate that the field of strain en¬ 
gineering in atomically thin semiconducting materi¬ 
als is very young as the current experimental works 
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FIG. 10: Lattice and electronic band structure of semiconducting TMDs. (a) Crystal structure of TMD 
where blue spheres represent a transition metal atom and yellow spheres, a chalcogen element, (b) DFT electronic band 
structure of single-layer M 0 S 2 and (c) bulk M 0 S 2 . Arrows indicate the gap in each case, (d) Two dimensional BZ. 


have achieved moderate strains in comparison to the 
maximum strains that 2D semiconductors can stand 
before rupture. At those larger values of strain, one 
would expect a strain-induced tunability of the band 
gap about 10 times larger than the maximum tunabil¬ 
ity experimentally observed so far. 


III. THEORETICAL STUDY OF STRAIN IN 
2 D CRYSTALS 

In this section we review the effect of strain on 
semiconducting 2D crystals from a theoretical point of 
view. We first describe the main features of the crystal 
and electronic structure of the most commonly studied 
2D crystals in the absence of external strain. Then we 
discuss the effect of homogeneous and inhomogeneous 
strain on the optical and electronic properties. As we 
have seen in Sec. |n| most of the current experimental 
results on strain engineering are performed on transi¬ 
tion metal dichalcogenides, and for strain levels that 
do not exceed ^ 3%. Here we include a revision of 
theoretical works on the field that make predictions 
for higher strain values than currently achieved ex¬ 
perimentally, as well as for other families of 2D semi¬ 
conducting crystals, like single and multi-layer black 
phosphorus, silicene and germanene, unexplored ex¬ 
perimentally so far. 


A. Electronic structure 

1. Transition metal dichalcogenides 

The most studied two dimensional crystals differ¬ 
ent from graphene are the semiconducting transition 
metal dichalcogenides, MX2 (where M =Mo, W and 
A = S, Se, Te). The crystal structure of monolayer 
and multilayer MX 2 is schematically shown in Fig. 
[lQ]^a). MX2 is composed of an inner layer of metal M 
atoms ordered in a triangular lattice sandwiched be¬ 
tween two layers of chalcogen X atoms lying on the tri¬ 
angular net of alternating hollow sites |lQQI41Q2l I1Q5F 
llQ7j . In Table |III| we give the values of the MX 2 lat¬ 
tice parameters. Similarly to graphene, the in-plane 
Brillouin zone (BZ) is an hexagon characterized by the 
high-symmetry points F = (0, 0), K= 47r/3a(l, 0), and 
M= 47r/3a(0, v^/2), where a is the in-plane M — M 
or X — X distance. 

The orbital character of the electronic bands of 
TMDs is studied in detail in Ref. [ml The conduc¬ 
tion and valence bands are made by hybridization of 
the dxy^ dx‘ 2 -y 2 and d^z^_x 2 orbitals of M, and the 
Px^ Py and Pz orbitals of X. All the single-layer MX 2 
compounds are direct gap semiconductors, with the 
gap placed at the two inequivalent K and K’ points 
of the BZ, as seen in Fig. [lQ|^b) for M 0 S 2 . The most 
relevant contribution at the edge of the valence band 
is due to a mixture of dxy and dx 2 -y 2 orbitals of the 
metal M, hybridizing with the Px and Py orbitals of 
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FIG. 11: Lattice and electronic band structure of black phosphorus, (a) Crystal structure of BP where blue 
spheres represent P atoms, (b) Top view of single-layer phosphorene. (c) DFT (dashed lines) and GW (solid lines) 
electronic band structure of phosphorene. The rectangular BZ is also shown. Adapted from Ref. cm with permission 
of the American Physical Society. 



a 

u 

c' 

M 0 S 2 

3.160 

1.586 

6.140 

WS 2 

3.153 

1.571 

6.160 

MoSe2 

3.288 

1.664 

6.451 

WSe2 

3.260 

1.657 

6.422 


TABLE III: Lattice parameters (in A units) of TMDs, as 
given in Ref. |10311105j . a is the M — M atomic distance, 
c is the distance between the M layers, and u the internal 
vertical distance between the M and the X planes. For 
bulk compounds, the z-axis crystal parameter is given by 
c = 2c'. 


the X chalcogen atoms. The main contribution to the 
edge of the conduction band is due to (i 3 ^ 2_^2 of the 
metal M, plus a minor contribution of Px an d Py or¬ 
bitals of the chalcogen X cm- In Table [lV| we give 


the energy gap values of the different compounds dis¬ 
cussed in this review. Multi-layers and bulk TMD 
compounds are indirect gap semiconductors, with the 
maximum of the valence band at F point and the min¬ 
imum of the conduction band at the Q point, which is 
situated midway between F and K points of the BZ. 


One of the main characteristics of TMDs is the 
strong spin-orbit coupling (SOC), which leads to a 
splitting of the valence and conduction bands. The 
valence band of M 0 S 2 and MoSe 2 is split by ^150 
meV, whereas the splitting for WS 2 and WSe 2 is con¬ 
siderably larger, of the order of ^ 400 meV, due to 
the heavier mass of W. The conduction band present 
a SOC splitting of several tens of meV at the band 
edge at the K point |102j . as well as at the Q point 


Band Gap 



Monolayer 

Bilayer 

Bulk 

Ref. 

M 0 S 2 

1.715 

1.710-1.198 

1.679 - 0.788 


WS 2 

1.659 

1.658 - 1.338 

1.636-0.917 

[T061 

MoSe2 

1.413 

1.424- 1.194 

1.393 - 0.852 

|106| 

WSe2 

1.444 

1.442 - 1.299 

1.407-0.910 

|106| 

ReSe2 

1.43 

— 

1.35 

|110| 

BP 

1.60 

1.01 

0.10 

|108| 

Silicene 

1.6 X 10“^ 

- 

- 

11091 

Germanene 

24 X 10-^ 

- 

- 

|109| 


TABLE IV: Band gap of 2D crystals, expressed in eV, 
as obtained from first principle calculations. Notice that 
all monolayer samples for M 0 S 2 , WS 2 , MoSe 2 and WSe 2 
present a direct gap. The two values given for bilayer and 
bulk samples for those compounds correspond to the sizes 
of the direct/indirect gaps, respectively. The gap of ReSe 2 
is direct for both, bulk and single layer. Results for BP are 
taken from GW calculations [108] . Results for silicene and 
germanene are taken from relativistic DET calculations, 
and corresponds to the low-buckled configuration (see text 
and Ref. ITO^ . 


ma, as shown in Fig. Elb). The phonon spectrum 
of TMDs is also very rich, and we refer the reader to 
Ref. |113j for a comprehensive review on the subject. 


2. Black Phosphorus 

Bulk black phosphorus (BP) is another layered van 
der Waals crystal consisting of puckered atomic sheets 
of phosphorus weakly coupled together |115j . The 
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Silicene 


Germanene 



FIG. 12 : Lattice and electronic band structure of silicene and germanene. (a) Side view of lattice geometry 
of low-buckled silicene and germanene. (b) Top view of the crystal structure, (c) First Brillouin zone of silicene and 
germanene. (d) Angle that defines the Si-Si bond and the direction normal to the layer, (e) and (f) show the relativistic 
band structures of silicene for 0 = 90° and 0 = 102°, respectively. The insets show a zoom near the K point of the BZ 
of the energy dispersion, with the gap induced by SOC. (g) and (f) are the same than (e) and (f), but for germanene. 
Adapted from Ref. ESI with permission of the American Physical Society. 


electronic and elastic properties of BP have been also 
studied in its nanotube structure ms]. A single layer 
of BP consists of a corrugated arrangement of P atoms 
with a thickness of ^ 5 A [see Fig. ll[a)-(b)]. Bulk 


BP is obtained by alternate stacking of the layers 
along the [001] direction. The BP intra-layer bond¬ 
ing is due to the sp^ hybridization of P atoms. 

The band structure of single-layer BP is shown in 
Fig. [^c) as calculated from DFT and and GW meth¬ 
ods. The orbital contribution to the valence and con¬ 
duction bands close to the gap has been studied in 
Ref. |108j . For the valence band, the orbital com¬ 
position at the F point can be written as |4^^(F)) = 
0.17\s) -\-0A0\px) -\-0.90\pz) whereas for the conduction 
band we have |4^‘^(F)) = 0.57|s) + 0.44|pa,) + 0.69\pz). 
The Pz orbital has the highest contribution in both 
bands, but the s and Px characters are not negligible. 
Interestingly, the Py orbital does not contribute to the 
formation of the valence and conduction bands. 

The gap at the F point obtained from GW cal¬ 
culations is ^ 1.60 eV for single-layer BP |lQ8j . 
which agrees well with photoluminescence measure¬ 
ments uni. As we increases the number of layers, 
the gap decreases monotonically reaching the value of 
^ 0.1 eV for the bulk material |108j (see Table IV). 
The in-plane electronic and optical properties of black 
phosphorus are highly anisotropic, due to the charac¬ 


teristic puckered structure of this crystal fTCTITT^. 


3. Silicene and Germanene 

Silicene is another 2D crystal that is attracting high 
interest. Encapsulated silicene field effect transistor 
have been fabricated recently, showing Dirac-like am- 
bipolar transport m- This material consists of a 
honeycomb lattice of Si atoms, similar to graphene, 
as shown in Fig. However, the two triangular 

sublattices made by the Si atoms are vertically dis¬ 
placed by 0.46 A, forming a buckled structure that 
resembles that of phosphorene. The interatomic dis¬ 
tance of silicene is 2.28 A, larger than that of graphene 
due to the large ionic radius of Si as compared to C. 
One can characterize the buckling of the 2D lattice 
by means of the angle 0 between the Si-Si bond and 
the direction normal to the plane, as sketched in Fig. 
p!^ d). Different configurations can be defined, like the 
sp^ (planar) for 0 = 90°, the low-buckled which corre¬ 
sponds to ^ = 101.73°, and the sp^ configuration, that 
corresponds to 0 = 109.47° nng. Due to the buckling 
of the lattice, the carriers move in a hybridization of 
the Pz orbitals with the a orbitals, leading to a large 
SOC as compared to graphene. As a consequence, an 
appreciable energy gap can be opened at the Dirac 
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points (see Table IV). 

Although it has not yet synthesized in experiments, 
germanene^ a single layer of Ge atoms, is another 2D 
crystal that it is attracting a lot of interest. The crys¬ 
tal structure is similar to that of silicene, with Ge-Ge 
distance of 2.42 A. The estimated gap, based on rel¬ 
ativistic DFT calculations, is also higher than in sil¬ 
icene (see Fig. [T^g)-(h)). In both cases, hydrostatic 
strain is expected to increase the size of the gap, what 
can help for the observation of quantum spin Hall ef¬ 
fect in the low-buckled structure of these compounds 

nog. 


B. Homogeneous strain 


1. Transition metal dichalcogenides 


culations for this kind of strain result in a lowering of 
the size of the band gap, as in the uniform case, with 
a possible transition to a metallic state. The electron 
mass at the Q point of the conduction band is found to 
increase (decrease) with uniaxial tensile (compressive) 
strain, while the effective masses at K decrease (in¬ 
crease) with increasing uniaxial tensile (compressive) 
strain um- The effect of mechanical deformations on 
the electronic band structure is especially relevant at 
the K point of the BZ, where both edges of the con¬ 
duction and valence bands are significantly lowered as 
compared to the unstrained case [53l II35j . On the 
other hand, the top of the valence band at F moves 
upwards. This is due to the fact that the orbital over¬ 
lap between the metal and chalcogen atoms decreases 
due to the change in the bond distances with applied 
strain [53] . 


The effect of homogeneous strain in the band struc¬ 
ture of 2D semiconducting crystals has been the focus 
of many theoretical papers. For TMDs, a number 
of DFT calculations have studied the dependence of 
the band structure and phonon modes with external 
strain [SU EH [Ml IMHISni, and the consequences 
on different physical properties have been discussed 
in Refs. |I31I[I33j . In Fig. it is shown the evo¬ 
lution of the band gap with uniform isotropic strain 
for single layer M 0 S 2 (a) and WS 2 (b). In the ab¬ 
sence of strain, TMDs are direct gap semiconductors, 
as discussed in Sec. Ill A I[ with the gap placed at 
the K point of the BZ. Uniform (uniaxial or biaxial) 
tensile strain leads to a linear decrease of the band 
gap, and a direct-to-indirect gap transition is reached, 
with the valence band edge being now situated at the 
F point [134] . Further tensile strain eventually closes 
the gap, leading to a semiconducting-to-metal transi¬ 
tion for elongations of the order of 11% [53l 11211 Il22j . 
This transition is expected to be possible in the lab, 
since for such strain the metal-chalcogen bonds are 
stretched, but not yet broken [53]. The case of com¬ 
pressive strain is similar, although its influence in the 
band structure is smaller as compared to tensile strain, 
as it can be seen in the top panel of Fig. [T^d). Notice 
that, due to the relatively small bending rigidity, com¬ 
pressive strain is difficult to be applied to suspended 
TMDs samples, although substrate-induced compres¬ 
sion is possible by depositing layers of TMDs on an 
elastomeric substrate, as discussed in Sec. |TTb 1 
Application of uniaxial tensile strain, on the other 
hand, is possible in the laboratory by the use of dif¬ 
ferent techniques, as reviewed in Sec. jir^ DFT cal- 


2. Black Phosphorus 

The effect of strain on the electronic and mechanical 
properties of single layer black phosphorus have been 
studied using first-principle methods nsa nsi nia- 
1146] . A detailed comparison of the BP electronic band 
structure using different numerical methods has been 
presented in Ref. 11361 Inclusion of many-body ef¬ 
fects within GoIFo increases the band gap as obtained 
from GGA-PBE. It was found that the electronic and 
the optical gaps, as well as the exciton binding energy 
increases (decreases) with the application of tensile 
(compressive) strain |I36[ 1145] . This trend is shown 
in Fig. T^a). The effect of uniaxial strain is simi¬ 
lar, as shown in Fig. [I4|^b), with a direct-to-indirect 
gap transition predicted for compressive strain of the 
order of ^ 5% [137] . Vertical strain applied per¬ 
pendicular to the BP layer has been studied in Ref. 
11381 This kind of strain leads first to a direct-to- 
indirect band gap transition, and for larger strains 
also a semiconducting-to-metal transition is predicted 
for uniaxial stress of the order of 24 GPa. 

BP has been proposed as an good material for ther¬ 
moelectric applications. Usually the dimensionless fig¬ 
ure of merit ZT = aTS‘^l{K.e + i^p) quantifies the per¬ 
formance of a thermoelectric material, where a is the 
electrical conductivity, T is the temperature, S is the 
Seebeck coefficient, /^e is the lattice thermal conduc¬ 
tivity and Up is the electronic thermal conductivity. 
Grystals with a large Seebeck coefficient and electri¬ 
cal conductivity and/or low thermal conductivity are 
expected to behave as good thermoelectric materials. 
Uniaxial strain has been proposed as a route to en- 
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Compressive strain 



FIG. 13: Homogeneous strain in TMDs. Evolution of the band gap of single layer M 0 S 2 (a) and WS 2 (b) under 
uniform isotropic tensile and compressive strain. The edges of the valence and conduction bands evolve with strain, and 
the points of the BZ at which they occur are indicated. The red circles mark the equilibrium structures in the absence 
of strain m- ( c) Dependence of the phonon frequency of the E^g and Aig Raman modes, as well as their difference, 
versus applied strain (tensile: red and magenta; compressive: green and blue; difference: yellow and cyan) [121] . (d) 
Evolution of the band gap with biaxial compressive (top) and tensile (bottom) strain for single-layer and bilayer M 0 S 2 
ED- Panels (a) and (b) adapted from Ref. [53] with permission of the American Physical Society. Panel (c) adapted 
from Ref. [121] with permission of Springer. Panel (d) adapted from Ref. |51j with permission of Elsevier. 


hance the Seebeck coefficient and the electrical con¬ 
ductivity in phosphorene, suggesting strain engineer¬ 
ing as an effective method to increase the thermoelec¬ 
tric performance of BP |139[ I141j . 

First principle calculations have found that biax¬ 
ial strain enhances significantly the electron-phonon 
coupling in black phosphorus. This effect has been 
proposed to favor superconductivity, leading to an 
increase of the superconducting critical temperature 
from 3 K to 16 K at typical carrier densities of 
^ 3.0 X 10^^ cm“^. This phenomenon is attributed 
to the simultaneous increase of the density of states 
around the Fermi level and the softening of the phonon 
modes. Furthermore, the anisotropy of the phospho¬ 
rene lattice leads to different effects on the electron- 
phonon coupling under uniaxial strain applied in dif¬ 
ferent directions, leading to a Tc ~ 10 K and 8 K for 
uniaxial strain applied along the x- and ^-directions, 


respectively 


3. Silicene and germanene 

The band structures of silicene and germanene are 
also affected by homogeneous strain pIM[T48l[T4^ . 
The results show that the size of the gap at Dirac 
points grows with compressive strain, and decreases 
with tensile strain, as shown in Fig. It is also 

observed that the gap is higher for higher angle 6>, as 
defined in Fig. Hid)- On the other hand, the Fermi 
velocity is only weakly affected by strain, being its 
value slightly less than the typical vp for graphene, 
this is ^ lO^m/s, due to the larger Si-Si interatomic 
distance. 

The characteristic buckled structure of silicene and 
germanene leads to a strain-induced self-doping phe- 
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FIG. 14: Homogeneous strain in BP. (a) Evolution of the electronic (Egap) and optical (Eopt) band gaps of single 
layer BP with biaxial strain. The dependency of the excitonic binding energy (Eexc) is also shown. The different 
numerical methods used are indicated as superscripts |136j . (b) Band gap evolution with uniaxial strain |137j . (c) Band 
gap with vertical compressive strain. The thickness of the unstrained layer is 2/io |138j . (d)-(e) Dimensionless figure of 
merit ZT as a function of carrier density for different values of uniaxial strain applied along the zigzag (d) and armchair 
(e) directions. Panel (a) adapted from Ref. [136] . panel (c) adapted from Ref. [138] and panels (d) and (e) adapted 
from Ref. [139] . with permission of the American Physical Society. Panel (b) adapted with permission from Ref. [137] . 
copyright 2014 American Chemical Society. 


nomenon, as predicted by first-principle calculations 
of Ref. [150] . The Dirac point is found to move be¬ 
low the Fermi level for compressive strain, leading 
to n—doped samples, whereas p—type doping can be 
achieved under tensile strain. 


C. Inhomogeneous strain 

In the previous section we have seen how homoge¬ 
neous strain can be used to control the electronics and 
optical properties of 2D semiconductors. In this sec¬ 
tion we will see how strain engineering can be used to 
modify and control those properties locally, by means 
of application of non-uniform strain. The possibility 
to continuously control the energy band gap in semi¬ 
conductors is of high interest for their application in 
optoelectronics devices. In particular, an inhomoge¬ 
neous strain field acting on a 2D semiconducting crys¬ 
tal leads to a continuous spatial variation of its elec¬ 
tronic band structure [94]. This effect has been pro¬ 
posed for the design of solar energy funnel, in which 
a broad band of the spectrum can be absorbed. As 
discussed in Sec. jll C 3[ this possibility was proposed 
theoretically by Feng et al. [m] and demonstrated 


experimentally by Castellanos-Gomez et al. in Ref. 
|94|for M 0 S 2 . 


A different route to achieve an exciton funnel has 
been proposed in Ref. [152] . where it was shown, 
based on first principle calculations, that a small twist 
between two semiconducting atomic layers creates an 
internal stacking translation i4(r) that varies with po¬ 
sition r and controls the local electronic bandgap 
Eg{u(y)). The authors of Ref. [152] propose the use 
of carbon/boron nitride or black phosphorus bilayer, 
and predict an energy gap with multiple local min¬ 
ima (which can act as funnel centers). The possibil¬ 
ity to locally tune the band gap in non-planar single 
layer black phosphorus has been studied theoretically 
in Ref. [153] . In that paper it was found that the gap 
can be reduced by 20% for non-planar phosphorene al- 
lotropes, with respect to its value on planar crystals. 
Triaxial in-plane strain has been proposed to engineer 
the gap in h-BN monolayer [154] . 
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FIG. 15 : Homogeneous strain in silicene. Depen¬ 
dence of the electronic energy gap (black circles) and Fermi 
velocity vf (red diamonds) near the Dirac points calcu¬ 
lated from relativistic first-principles method, under dif¬ 
ferent applied strain. The inset shows the energy of unit 
cell for different hydrostatic strain. Adapted from Ref. 
[ 109 ] with permission of the American Physical Society. 


1. Tight-binding approach considering inhomogeneous 
strain: application to TMDs 

Theoretically, the study of non-uniform strain by 
means of first-principle methods is extremely expen¬ 
sive in computational terms, due to the huge unit cell 
that must be considered. An affordable alternative is 
the use of tight-binding methods |155j . which are a 
good base to further include electron many-body ef¬ 
fects, electron-lattice interaction, etc. Using geomet¬ 
ric arguments, the energy hopping integral Vpg be¬ 
tween two orbitals can be written in terms of Slater- 
Koster two-center parameters Vssa^ Vspa^ Vppa^ ^ppir^ 
Vddcj, VddTT, Vdd6, Vpda, VpdTT^ etc., which, together with 
the crystal fields, represent the suitable tight-binding 
parameters to be optimized m- 

Within the Slater-Koster framework, hopping pro¬ 
cesses between two atoms depend only on the relative 
angle and on the relative atomic separation Rijpi, = 
\Rij^u\, through the hopping integrals 
where i,jf label the unit cell and the /i, indices de¬ 
note the specific atom and orbital within the unit 
cell. When considering the effects of homogeneous 
and inhomogeneous strain it is common to assume, 
following Refs. [15611157] . that the main contribution 
comes from the modulation of the hopping integrals 
distance, neglecting the weak 
dependence on the relative angles. 


For a given profile of the strain tensor e{x^y^z)^ 
assuming the characteristic length of the sample to be 
much larger than the lattice constant, we can write 


'R^ijpiv{p^ 1 V-) 


I-hi(x,y) 


p>0 


( 2 ) 


where is the relative interatomic distance in the 

absence of strain. Therefore the hopping integrals will 
be modified as 





1 Pi,j,fi,iy 

R^ 

y) 

(3) 

- 





jdO 

^ijpiV 

The 

parameter 



= 


[Rijy^] 


where 5Rijpy{x,y) = 
deformation coupling 
—d\ogVi^j^p^j^{R)/d\ogR characterizes the electron- 
phonon interaction strength corresponding to each 
hopping bond, which can be evaluated microscopically 
from first principle techniques. 

It is illustrative to consider the modification of the 
band structure of a M 0 S 2 ribbon with a non-uniform 
profile of strain, originated from delamination from 
the substrate [94]. For that aim, it is useful to con¬ 
sider a one-dimensional profile along the ^-axis of the 
strain tensor. Neglecting, as a first approximation, the 
change in the membrane thickness with the in-plane 
applied strain, one can write a strain tensor 


£(.x,y,z) = i{y)=en 



(4) 


where L is the characteristic width of the ribbon, 5max 
the maximum strain attained at the center of the rib¬ 
bon y = L/2, and z> is the Poisson’s ratio of the mem¬ 
brane and f{y) is the function that accounts for the 
distribution of strain along the sample. For a half- 
sinusoidal distribution we have f{y) = sin(^), but 
other distributions of the strain along the sample, like 
Lorentzian distribution, Gaussian distribution, etc., 
can be considered with this model 11581. The results 


are shown in Fig. 16 Panel (a) shows the band struc¬ 
ture obtained for different levels of strain. One ob¬ 
serves that the direct band gap is reduced by appli¬ 
cation of strain, in agreement with experimental ob¬ 
servations, as discussed in Sec. |II C| The local change 
in the gap due to non-uniform strain is shown in a 
remarkable way by the bending of the local density 
of states (LDOS) along the wrinkle, as shown in Fig. 
[T^b). The plot shows how the gap is reduced at the 
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FIG. 16: Inhomogeneous strain in TMDs. (a) Band structure for a M 0 S 2 ribbon with zigzag edges. Periodic 
boundary conditions have been used, and a non-uniform strain profile as given by Eq. Q, where £max correspond to 
0%, 2% and 4%, respectively. The arrows indicate the direct transitions at the K point of the BZ, which are observed 
by photoluminescence experiments, (b) Color map of the local density of states of a 200 nm M 0 S 2 wrinkle, as a function 
of the position along the wrinkle. The strain profile is sketched in the top panel, (c) LDOS calculated at the center 
X = 100 nm (maximum strain) and at the edge x = 0 (unstrained). The arrows denote the estimated band edges as given 
by the second derivative method. Adapted with permission from Ref. 1941 copyright 2013 American Chemical Society. 


center of the ribbon, for which the strain is maximum, 
as it is seen by the cuts shown in Fig. [T^c). 

2. Gauge fields and topologieal aspeets 

Local modulation of the strain can be used to mimic 
gauge fields in the samples |159j . This opens the 
way for studying more fundamental phenomena in 2D 
crystals, for example, realizing giant pseudo-magnetic 
fields and associated zero-field quantum Hall effects 
as observed in graphene 

Cazalilla et al. have suggested in Ref. 11611 that 
time-reversal invariant topological phases and quan¬ 
tum spin Hall effect can be induced in single-layer 
TMDs as M 0 S 2 (at low carrier densities) due to the 
lack of inversion symmetry and SOC. Indeed, semi¬ 
conducting TMDs under shear strain are expected to 
develop spin-polarized Landau levels in different val¬ 
leys. For M 0 S 2 , gaps between consecutive Landau lev¬ 
els have been estimated to scale as hujc/ks ^ 2.7Ho[T] 
K, in terms of the cyclotron frequency ujc = ^ 

where m* is the effective mass, and Ho[T] is the 
strength of the pseudomagnetic field in units of Tesla. 
Taking into account that, for graphene, pseudomag¬ 
netic fields of the order of ^ 100 T have been already 
observed experimentally and considering the in¬ 

trinsic limitations imposed by the maximum tensile 
strength and the value of the shear modulus of these 
materials, gaps between Landau level of up to ^ 20 


K are within experimental reach. A different route to 
tune the electronic topological properties is the appli¬ 
cation of an electric field, as discussed in |162j . For a 
comprehensive review of effects of strain in graphene 
and other 2D crystals, from a more fundamental point 
of view, we refer the reader to B. Amorim et al. una. 


IV. CONCLUSIONS AND OUTLOOK 

We reviewed the effect of strain on the new fami¬ 
lies of semiconducting 2D crystals. The main parts 
of the review are devoted to the techniques to fab¬ 
ricate and induce strain in these materials, and how 
strain engineering can be used to modify in a con¬ 
trolled manner their optical and electronics proper¬ 
ties. Especial attention has been paid to transition 
metal dichalcogenides (M 0 S 2 , WS 2 , etc.), since the 
reported experiments are only limited to this family 
of materials so far. Several proposal to use band gap 
engineering for optoelectronics and photovoltaics ap¬ 
plications have been discussed. In particular, the use 
of localized strain can make it possible to fabricate so¬ 
lar energy funnels which can be used for photovoltaic 
cells that absorb a broad band of the solar spectrum. 

Apart from the modification of the electronic band 
structure, straintronics can exploit the piezoelectric 
properties of some families of 2D crystals, converting 
mechanical to electrical energy. Piezoelectricity is a 
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well-studied effect in which stretching or compressing 
a material causes it to generate an electrical voltage, 
or the reverse, in which an applied voltage leads to 
contraction or expansion of the material. M 0 S 2 , a 
material which is not piezoelectric in its bulk configu¬ 
ration, has been shown to become piezoelectric when 
it is thinned down to a single atomic layer (or to an 
odd number of them), and flexing it in the proper di¬ 
rection, as it was shown in Ref. |164j . Interestingly, 
the authors monitored the conversion of mechanical 
to electrical energy, and observed voltage and current 
outputs, finding that the output voltage changed sign 
when the direction of applied strain is reversed, as 
expected due to lack of inversion symmetry in single 
layer TMDs. 

Since 2D crystals can be stretched much farther 
than conventional materials, these observations can 
open the door to development of new applications 
for these materials and its unique properties. It is 
also worth mentioning that, although in this review 
we focused on the tunability of electronic and optical 
properties, strain-induced effects can be employed also 
to tune alternative degree of freedoms, as magnetism 
|165l I166j or lattice dvnamics ns?]. 
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